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A recent experiment with squeezed hght has shown that two-photon absorption by an 
atom can occur with a linear intensity dependence. We point out that this result verifies 
a prediction made by us a decade back from an analysis of a non-local model of Quantum 



^ \ Electrodynamics. This model had earlier been proposed by us in an ad hoc manner to 

O ■ interpret certain features of multiphoton double ionisation and above-threshold ionisation 

Q . in an atom placed in a strong laser field ; in this paper we show that the model can be 

0> ! obtained field-theoretically by demanding covariance of the field Lagrangian under a non- 

- ' local U(l) gauge transformation. We also show that the model makes direct contact with 

squeezed light, and thus allows us to describe these two completely different scenarios from 
p ■ a unified point of view. We obtain a fundamentally new result from our non-local model 

1^ , of QED, namely that only the past, hut not the future, can influence the present — thus 

establishing a non-thermodynamic arrow of time. We also show that correlations within 
a quantum system should necessarily he of the Einstein- Podolsky- Rosen (EPKj-type, a 
/S ' result that agrees with Bell's theorem. These results follow from the simple requirement 

c^ ■ of energy conservation in matter-radiation interaction. Furthermore, we also predict new 

and experimentally verifiable results on the basis of our model QED. 
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1. INTRODUCTION 

In a beautiful experiment, Georgiades et al [1] generated a squeezed vacuum field by 
an optical parametric amplifier and fo cussed it onto a cloud of cesium atoms in a magneto- 
optic trap. They found that the rate of the two-photon transition 6S'i/2 -^ 6D5/2 deviated 
from the expected quadratic intensity dependence to a form R = al + /3/^, where / was 
the squeezed vacuum intensity ; the linear term dominated at low intensities, and the 
quadratic term at higher intensities. Our conventional wisdom tells us that this feature 
of linear intensity dependence of two-photon transitions is shared neither by thermal light 
nor by coherent light. 

In this paper we draw the reader's attention to a non-local model of Quantum Elec- 
trodynamics (QED), suggested quite a few years ago by us [2] in which a possibility of this 
linear intensity dependence of two-photon transitions was pointed out for the first time, 
albeit in a different context of multiphoton transitions and above-threshold ionisations in 
an atom placed in a laser field. We were led to this model of non-local QED through an 
analysis of the behaviour of an atom placed in an intense laser field under certain specified 
conditions — namely, that two distinct photon absorption events by the atomic electron/s 
be correlated in phase — and we obtained this two-photon transition with a linear intensity 
dependence as the central result of the paper. In [2], however, the non-local potential was 
introduced in an ad hoc manner ; in the present paper we provide a field-theoretical basis 
for the same and briefiy repeat the analysis made in [2] . We bring out certain parallelisms 
between the two contexts (intense laser field and squeezed vacuum field) and show that 
the aforesaid model of non-local QED can be applied to both. As a particular result of 
this application, we find that the time evolution of an (atom -|- field) system under the 
said conditions shows a remarkable resemblance to a squeezed vacuum. (As we shall show, 
apart from pre-multiplying factors, they turn out to be identical). We further show that 
a specific result for a single-electron atom placed in a non-local electromagnetic field ex- 
actly matches the abovementioned experimental result of two-photon transition in Cesium 
placed in a squeezed vacuum. 

Subsequent to our work [2], a few studies on squeezed light (e.g. Gea-Banacloche 
[3], Javanainen and Gould [4]) predicted this linear intensity dependence from certain 
special models of the four-field correlation function in squeezed light fields. Our approach 
needs no such squeezed-field-specific modelling, but is valid for all (classical as well as 
non-classical) radiation fields where phase-correlated photon absorption/emission can take 
place in interaction with matter, as we shall show below. 

Consider an intense classical (thermal or laser) light incident on an atom, and let a 
pair of photon absorption events occur at two different space-time points (r, t) and (r', t'), 
the interval between which is very close to zero, i.e. |r — r'| ~ c|t — 1'|. This situation can 
arise when, say, the outer electrons of a Rydberg atom absorb a pair of photons within a 
time 5t <^ —. For laser intensities of order lO^'^ W/cm^ or higher, this is commonplace, 
if we remember that in the visible range, a fiux of one photon Uq tQ corresponds to ~ 
5 X 10^^ W/cm^ of intensity (oq and to being the atomic units of length and of time 
respectively). We have here a situation where two events, which may have occurred far 
apart in space, are correlated in phase ; indeed, if the events are spatially separated, the 



correlation assumes the character of an EPR-hke "quantum entanglement" . 

Consider now a different scenario : radiation from an optical parametric oscillator 
generating a phase- squeezed vacuum field is absorbed by a three-level atom having a prin- 
cipal quantum number > 5. If the squeezed vacuum is produced by parametric frequency 
down-conversion, where an incident pump photon has split up into two lower-frequency 
photons (signal and idler) which constitute a highly correlated photon pair, a two-photon 
absorption similar as observed in [1] occurs ; i.e. two photons in exact (or nearly exact) 
phase correlation can be absorbed in two events widely separated in space even at low 
intensities. 

The feature common to these two entirely different scenarios is evident. The latter 
scenario is directly linked to the abovementioned experimental finding in [1], while the 
former is linked to phase-correlated multiphoton excitation/ionization in a regime where, 
however, perturbation theory holds. In this paper we show that our model in [2] can be 
used, by exploiting this common feature of correlated atomic response to radiation at two 
points far apart in space, to give a theoretical description of these two different scenarios 
from a unified viewpoint. As already mentioned, the linear intensity dependence of two- 
photon processes comes about as a natural result of this model, and we shall show how this 
linear intensity dependence serves to explain several features long known in multiphoton 
double ionisation as well as in above-threshold ionisation, and also serves to predict some 
new results. 

2. THEORY 

In the usual simple harmonic oscillator mode expansion of the electromagnetic field 
potential A, an integral number of photons (0,1,2,...) is associated with each mode. All 
these modes are orthogonal and independent of each other, and since photons obey the 
Bose statistics, each mode can be occupied by any arbitrary number of photons. However, 
in situations such as described above where a phase- correlated behaviour of photons be- 
comes manifest, one might think whether a different, correlated type of mode expansion 
could better describe the situation, much as a correlated wave-function gives a better rep- 
resentation of, say, the Helium atom than the orbital product wavefunction does. Since 
there is no force acting between photons, the observed "correlation" between them is not 
due to any "configuration interaction" as in the case of a multi-electron system, but due 
either to the presence of certain non-classical features in the radiation field introduced by 
a non-linear medium (e.g. squeezing), or to correlated interaction of the radiation field 
with the atomic electrons, or both. 

We start from the basic fact that in both these scenarios the radiation field (either as 
generated, or as detected by the atomic electrons) enjoys a large second-order coherence, 
i.e. two photons are coherently absorbed by the atomic electrons within a vanishingly small 
phase difference Acj) = ivdt. In a study of two-photon transition using a phase-squeezed 
light within the framework of the usual QED, Javanainen and Gould [4] observed that 
the "four-field correlation function G'(2)(r, t; r', t') = {E-{r,t)E-{r\t')E+{r\t')E+{r,t)), 
essentially the joint probability that two detectors placed at r and r' will record photon 
counts at times t and t', does not always suffice to predict the two-photon transition rate. 



In fact, we need a more general correlation function. . . We resort to heuristic modelling of 
the field correlations" . 

In both [3] and [4] the authors gave new four-field correlation functions derived from 
specific properties of squeezed light. In the present paper we show that neither the use 
of such specific properties nor any heuristic modelling is necessary, if we adopt the earlier 
proposal of attempting a correlated description of the electromagnetic field. 

We show below that the non-local model of QED as expounded in [2] serves to bring 
out the feature of (phase) correlation between two (or more) photon modes. Consider 
two zero-order dressed states |\1//) = |'0/)|n/) and |\l/i) = \'^i)\ni). In the Coulomb 
gauge the radiation interaction matrix element coupling these two states is Mfi{t) = 
— ^(•0/(r),n/|p.A(r, t)|'0i(r), rii). As in our earlier work [2], we introduce a non-local 
potential ^(r, t, r', t') such that the above matrix element becomes 

g 
Mfi{t) = {i(jf{x),nf\p.A{x,x')\iJi{x'),ni) where x = {r,t), x' = {r',t'), (1) 

i.e. the interaction at the point x is correlated with the interactions at all other points 
x' . This non-local interaction had been defined in [2] in an ad hoc manner ; in this paper 
we show (see Appendix A) how such a potential can be introduced field-theoretically by 
demanding covariance under a non-local U(l) gauge transformation. 

The integration in (1) runs over dr and dx' , and in this paper we set the limits of the 
time integral dt' either (i) from — oo to t — ' (giving the retarded interaction), or (ii) 

from t + '^~^ ' to oo (which gives the advanced interaction), leaving out EPR "quantum 
entanglement" -like correlations for the present. We shall discuss the question of causality 
later and shall find that this restriction turns out to be unnecessary. 

We can equivalently use a similar non-local form of the multipolar interaction at the 
space-time point x : 

—e / Y.£{x,x') dx' {!') 

Although nothing can be said about the specific mathematical forms of these non-local 
potentials/field strengths, we will find that we can go a long way to predict specific conse- 
quences of this general non-local field. In order to quantise this non-local field, we proceed 
by carrying out a Fourier expansion [5] of this non-local potential over the usual photon 
modes : 






A{x, ^') = 2_^ 2^ [Ci2ei2(r, r) exp(zki.r + zk2.r' - iuj\^^t - iuj\^^t') + c.c] (2) 

kiAi k2A2 

where we have written C12 = C\^^\^\^^\^, and 612 = ekiAik2A2- ^^ view of (1), we cannot fix 
the direction of the polarisation vector e(f, r') "by hand" (as in standard QED) without 
loss of generality ; this is at once obvious if we make the dipole approximation in (1) or (1'). 
Later in this paper we shall use for e(f, f') a power series in f and f' with undetermined 
coefficients, which is the most general form possible. However, to quantise the field using 
the Coulomb gauge we impose for the present the condition that 

(ki + k2).e = 0. (3) 



This restriction will be removed later. 

From eqn. (A2) in Appendix A, we note that the non-local electric field in the Coulomb 
gauge is given by 

£(a:,xO = -(l/c)('|^ + ^') =^5] 5^K,+a;kJ[C'i26i2(f,f')/(x,x')-c.c.] (4) 

kiAi k2A2 

where f{x,x') is the phase factor in (2). The expression for the non-local magnetic field 
can be found in an analogous manner. 

To obtain the field energy at a point (r, t) we need E'^ and B^ which we define similarly 
as in (1) ; thus 

£;2(r,t)= f [ £{x,x')£{x,x")dx'dx" 

where the limits of integration over dt' and dt" are determined by the choice of interaction 
(retarded or advanced). However, in the present work r's are of atomic dimensions, and we 
simply put t', t" < t for the retarded interaction and t', t" > t for the advanced interaction. 

An examination of eqns. (2) and (4) shows that in E'^ and B^, the crossed terms be- 
tween two modes kiAi and k2A2 vanish unless the two modes are strongly phase-correlated. 
Denoting the population fractions of such correlated-pair modes and of the independent 
modes by an and aj respectively, use of standard methods in Quantum Electrodynamics 
leads to the following expression for the field energy in a volume O : 

W - ajWi + aj^Wii (5) 

where Wi is the usual energy summed over single modes, and 

^ii = iE(^k'A + 4c^^^L), (5') 

kA 



where YkA = -«^^-^2c^k(CkA - C^x)^ ^kA = ^^-^(CkA + C^^)- Each mode subscript 
here actually stands for a correlated mode pair (kiAi, k2A2), and 2a;k stands for cuki +'^k2- 
Y and Z satisfy the canonical equations of motion with Wu as the Hamiltonian as in 
standard QED, and we now quantise this non-local field by requiring that Y and Z be 
g-numbers obeying the commutation relation 

[^kA, YWX'] = ^^^kk'^AA'- (6) 

We thus obtain the non-local field energy in terms of a new pair of creation and 
annihilation operators : 

mi = i E ^^k (^^kA^A + ^A^kA) (5") 

kA ^ ^ 

where 6kA, 4a = /^^ — (2c^k^kA ± «>kA) • (7) 
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Equation (5") is identical in form with standard QED, but the creation and annihilation 
operators obey a new commutation relation : 



^kA,\/A' 



2'5kk"^AA', [bkX.W'X'] 



"kA'^k'A' 







(8) 



as may be directly verified from (7) and (6). The following, and only the following, results of 
operation of these new correlated-pair mode operators are consistent with the commutation 
rules (8) : 

^kAl^^kA) = v^^^kA|?^kA - 2), &Ll^kA) = V?^kA + 2|?ikA + 2), bl^b]^x\n]^\) = ni^\\ni,x) 

This shows that b and b' are two-photon annihilation/creation operators, and b'b is the 
usual number operator. Once again, by creation of two photons we mean the creation of 
one photon in each member of a correlated mode pair, and similarly for annihilation. We 
can now write down the potential and field strengths in terms of these new operators ; 
thus 



A{x,x') = y^^\l ^ e]^x{f,r) b]^xf{x,x') + b1^f*{x,x') 



£{x, x) = i 



kA 

E 

kA 



OcUk 



2TvhuJ]f^ 



ekA(r, f') b]^xf{x,x') 



blxf*ix,x') 



(9a) 



(96) 



The results obtained so far from our non-locality postulate may be summed up as 
follows : quantization of a radiation field with a high degree of second-order coherence 
can depart from the usual QED and can lead to photon modes (more exactly: correlated 
mode pairs) generated by two-photon creation and annihilation operators. Indeed, if we 
let a zero-order dressed state of an atom \i(^)\n) evolve in the presence of the non-local 
electron-radiation interaction (1) or (!'), we obtain a linear superposition of dressed states 
of the form 



oo 



\^)yc,{t)\2k) 



fc=0 



(10) 



for even n ; for odd n, \2k) is replaced by \2k + 1). The possibility of emission/absorption 
of two photons in the first order of interaction, i.e. linear in field intensity becomes evident 
at once as a consequence of this model. Also, the commutation relations (8) are gauge- 
independent, and hence from now on we can do away with the transversality restriction as 
mentioned after eqn.(3). 

3. RESULTS AND DISCUSSION 

Before we work out the possible consequences of this non-local field in a typical atomic 
transition, let us note its formal similarity with a squeezed state. We recall the definition 
of the squeezed state [6] : 

|a,C)=I^(«)^(C)|0), (11) 
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where the squeeze operator S{Q = exp[i^*a^ — ^({a')"^], and the displacement operator 

2 J. 

D{a) = exp (aa' —a*a) = e 2 e"" e~" °'. {a and C are arbitrary complex numbers.) 
Putting a = gives the squeezed vacuum, while putting C = gives the coherent state. 
An expression for the squeezed vacuum in terms of the number states has been given by 
HoUenhorst [7] ; putting ( = re*^, this becomes 

00 — 

^(()|0) = (coshr)2 E i^e'Hanhr)'" ^ \2n) (12) 

n=0 

Thus we see that the squeezed vacuum 5'(C)|0) has the same form as the expansion 
(10) above. The two-mode squeezed vacuum [8], where the squeezing operator is given by 
the more generalised form 

S{Q = exp[C(c^)at(c^)at(20 - c^) - C*(c^)a(c^)a(20 - c^)] 

provides another example where the doubly occupied modes in (12) would be replaced 
by singly occupied correlated mode pairs ; the signal and idler photons generated in a 
parametric down-conversion present themselves at once as forming such a correlated pair. 
We can therefore say that the squeezing operator S{Q is mathematically identical (except 
for coefficients/multiplying factors) to the evolution operator in a non-local field, and hence 
atoms interacting with a squeezed field and with an intense radiation field (for which a 
non-local description is valid) can be expected to show the same non-classical features. We 
find below that it is indeed so. 

We now work out two typical first-order transition matrix elements using this non- 
local description for a radiation field. Some of these results have appeared before [2b] , but 
are included here for the sake of completeness. First it must be mentioned that depending 
on the photon fiux and on the coherence properties of the laser field, one must in general 
describe the field as 

S = aL^L + ttN^N (13) 

where the subscripts L and N refer to the local (usual) and the non-local fields respectively 
{aL and a^ are identical with ai and an in eqn.(5)) ; the relative strength of the two 
coefficients would depend, essentially, upon the degree to which the response of the atomic 
electron/s to the field appears phase- correlated. 

We work within the dipole approximation, and formulate the typical matrix element 
of this two-photon absorption operator (I) for a single-electron atom, and (II) for a two- 
electron atom. 

(I) For simplicity we confine ourselves to a single radiation mode |nkA)- Define 

Ti = {'^l)lJmf{x)■:nl,\-'2\v£N{x,x')\'^l)l^^rn^{x'),nl,x) (14) 

where the /'s are the orbital angular momenta in the initial and final states. We use 
equation (9b) for Sn{x,x') and note that the time integration in (14) for the non-local 
interaction yields the required energy balance Ef = Ei + 2?iuj\<^ if and only if we use the 



retarded interaction ; use of the advanced interaction fails to satisfy energy conservation. 
(See Appendix B for the mathematical details.) This remarkable result shows that causality 
as we understand it — ability of the past, hut not of the future, to influence the present 
— follows as a necessary condition for energy conservation in our non-local picture of the 
electromagnetic field ; we have obtained a non-thermodynamic arrow of time having a 
purely quantum nature. 

Incidentally, energy conservation is satisfied even if we set the upper limit of the 
integration over dt' (eqn.(B3), Appendix B) to t exactly, i.e. put c = oo, which corresponds 
to an EPR-like "quantum entanglement". Thus we find, as had been mentioned in the 
beginning, that such situations need not be left out from our non-local picture, and that 
causality in a quantum system need not be restricted by the special theory of relativity. 
However, if we go further in a closer analysis of the integration over dt' , we find that we 
can make an even stronger statement. Consider the case where |r — r'| > ct, in which case 
the upper limit in the t'-integral in (B3) becomes negative (remember that both r and r' 
extend over all space), and the integral vanishes, leading to non-conservation of energy ! 
Thus in our non-local picture of the quantum world, attempts to combine causality with 
the special theory of relativity are not only unnecessary, but can also lead to results which 
do not agree with experience ; one is immediately led to think of Bell's theorem. Note 
that such an EPR-like correlation between x and x' is also indicated by the commutation 
relation (A3) in Appendix A, which gives the necessary condition for gauge invariance of 
the non-local electromagnetic field tensor. 

Finally we get (in atomic units) 



Ti = ^^(V.;(r)|r.e(f,f')IV'^(r')) (15) 

As mentioned earlier, we use for e(f, f') a power series in f and f' with undetermined 
coefficients ; the specific form chosen is 

oo 

e{f,r') = {f + f')\^cinPn{^-^') (16) 

n=0 

(Symmetry demands that the coefficients of f and of r' be identical.) Putting 'i/'i(r) = 

]\j.RiXnYi.^.{f) and similarly for tpf gives 

Ti = iN.Nf X l^n{Ii + J2) (17) 

where 7^ is a radial integral, and the X's are two angular integrals, one of which gives the 
selection rule // = U, while the other gives // = U, U ± 2. (See Appendix C for details.) 
Thus from (17) we obtain two-photon transitions of the type S ^ S, S ^ D. . . linear in 
intensity. We also note that the squeezed vacuum corresponds to the limit a -^ in eqn. 
(11) ; an increase in a, which results among other things in an increase in intensity, leads 
to a deviation of the squeezed field from the form (12) to a situation where the summation 



over mode indexes runs over all integers (instead of only even or only odd integers as in 
(10)), so that in the first order, an increase in intensity in squeezed light-matter interaction 
causes, progressively, single- photon transitions instead of the two-photon transitions as was 
the case in the squeezed vacuum. This result exactly matches the observation of Georgiades 
et al [1]. Other possible two-photon absorptions linear in intensity in a squeezed vacuum 
can be trivially predicted. Two possibilities immediately come to mind that should be 
experimentally feasible : (i) two-photon ionisation of a Rydberg atom, and (ii) two-photon 
electron detachment from a negative ion, both with a squeezed vacuum. Indeed, a whole 
field of non-destructive two-photon spectroscopy of atoms/molecules opens up. 

At this point we should remark that the formal similarity between a squeezed field and 
a laser field with strong second-order coherence can become misleading when we consider 
the intensity variation of the fields interacting with matter. As we shall soon see, the higher 
the intensity, the more phase- correlated does the classical (thermal or laser) radiation field 
"appear" to two detectors (e.g. two atomic electrons), or one detector at two different 
times, rendering this non-local description the more valid. On the other hand, as we just 
saw, increase in intensity in a squeezed vacuum takes us back, so to say, into the classical 
domain (the two-photon transition becomes non-linear with increase in intensity). 

(II) For a two-electron atom we require a matrix element of the type 

Tu = {'^ifi'ix,x'),ni^x -2\{r + r').£N{x,x')\tjJl^^l'{x,x'),n]^x) (18) 

As before, the time integration in (18) yields the energy balance S{Ef, Ei + 2huj]^), leaving 

T„=^^(V;y(r,r')|(r + r').e(f,f')|V^(r,r')) (19) 

To evaluate Tu we use a correlated wavefunction of the form [9] 

i^^R-^^Yl ^a*(^)*m(^; «, ^ r) (20) 



where R = Vr^ + r'^ and a = arctan {r' /r) are the hyperspherical coordinates. 

Substituting eqns. (20) and (16) in (19) we obtain (see Appendix D for details) 



T„ = i^^ I F}{R)RF,{R)dR j AfA.P^n/'^ '\cos 2a) {cos a] 
X (sina)'/+^i(cosa + sinQ;)P^"*'^*^(cos2Q;) cos^ a siv? a da 



If+h 




Tli 

OO 



^fXz^(^'^')] {l + Pi{f.r')J2^kPk{f.f')yff^^i,{r,f')dfdf' (21) 



k=0 



We obtain the selection rules for If and l'^ from the two angular integrals, one involv- 
ing '^i^cikPki'f'.f') and the other involving '^i^akPi{r.r')Pk{f.r'). They are tedious but 
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straightforward ; note that although the coefficients a^ are unknown, the relative strength 
of the two terms can be analytically obtained. For the first integral, both A/ and A/' are 
even for even k and odd for odd /c, whereas for the second integral A/ and A/' turn out to 
be odd for even k and even for odd k. Thus finally we see that parity is conserved, as in 
the single-electron case. 

For carrying out the integration over da one can use the result [10] 

1 

(1 - x)P{l + xY P^^'''^\x)P^'^\x) dx 
■1 

2P+-+1(1 + a)„(l + 7)^r((7 + l)r(p + 1) y. {-m)r{l + 7 + ^ + Hr(l + P)r 



m\n\V{p + a + 2) ^ r!(l + 7)^(/9 + a + 2)^ 

X 3i^2(-n, l + a + /3 + ?i, l + p + r;l + a,2 + p + (T + r;l) (22) 

On symmetry grounds. An = Uf — rii is always even ; sample calculations show that the 
integral over da is strongly peaked when n/ = n^. We see therefore that this two-photon 
transition maintains the radial correlation in the two-electron atom. 

We can now summarise the properties of the matrix element T as follows : 

(i) It is a two-photon excitation linear in the field strength ; 

(ii) It conserves the parity of the atom ; 

(iii) For a two-electron atom, it conserves the degree of radial correlation. 

In view of eqn.(13) this double excitation will be running, in general, parallel to 
usual single-photon transitions in an experiment, and hence in any multiphoton excitation 
process in a field having a strong second-order coherence the overall transition would have 
the form TTDDTD ... or some such random mixture (D's are single-photon transition 
moments) , tending towards the extreme form TTTT ... as the field grows more intense. 
If we start out with a closed-shell atom in an S'-state, such a series of T-transitions would 
result in the formation of strongly correlated doubly excited states which are known to 
favour simultaneous double electron escape (Wannier [11] ; see also the review articles of 
Fano [9] and Rau [12]). This is to be contrasted with sequential two-electron escape, which 
does not involve any electron correlation. 

We now demonstrate briefiy that many of the features observed in multiphoton double 
ionisation in atoms placed in strong laser fields can be shown to follow as natural conse- 
quences of this non-local field model. A fuller demonstration can be found in [2b]. As 
mentioned at the outset, interaction of matter with a strong radiation field fits into this 
non-local picture if two photons are absorbed by an atom within a time interval 5t <^ —. 
For practical purposes this requirement can be replaced by a more convenient one, namely 
that in a time — , the number of photons fiowing into and out of the "volume" V of the 
atomic shell occupied by the outer electron/s be much larger than the number of these 
outer electrons. Thus for an n-(outer)-electron atom the condition reads 

^ > n 23 
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On the other hand, when the correlated motion of two electrons is involved (see discussion 
after eqn.(22)), we have to remember that the laser electric field should be significantly 
smaller than the Coulomb fields within the atom, in order that the electrons can maintain 
their mutual correlation. Choosing a field of order 10~^ a.u. gives us a crude estimate for 
an upper limit 

/max ~ 10^' W/cm^ (24) 

to the laser intensity upto which we can expect electron correlation in a two-electron atom 
to be maintained. 

We now examine some physical consequences of our model QED and compare them 
with observed features in Multiphoton lonisation (MPI) and Above-threshold lonisation 
(ATI) experiments. Earlier, the occurrence of a TT. .-like chain was shown to lead to the 
possibility of direct double ionisation (as opposed to sequential double lonisation) [13] ; 
we now see that the possibility can materialise if and only if both the conditions (23) and 
(24) are fulfilled. Simple considerations show that for short wavelengths the two conditions 
can become mutually exclusive, thus ruling out the possibility of Direct Double Ionization 
(DDI), whereas for long wavelengths both conditions (23) and (24) can be satisfied over 
a broad range of intensities. There is a caveat, though ; for very large wavelengths the 
condition (23) can be satisfied at quite low intensities, but the rather large number of D's 
in the multiphoton excitation chain TTDDTD . . . will tend to spoil the electron correlation 
pattern. Also, for a given power output, the peak intensity in short laser pulses is higher 
than in longer ones, and hence we should expect DDI to occur more readily with shorter 
pulses than with long pulses. This accords with the observations of L'Huillier et al [14] 
who measured multiphoton ionisation in Xenon using both nanosecond and picosecond 
pulsed lasers at 1064 nm and 532 nm, and those of Agostini and Petite [15] and Delone 
et al [16] who measured MPI of Ca at 1064 nm with picosecond pulses and nanosecond 
pulses respectively. Furthermore, Yergeau et al [17] found that in MPI of rare gases using 
a CO2 laser (A = 9.55/U and 10.55/u), the multiply charged ions are formed sequentially. 

In a strong laser field, another interesting effect stems from the photon number-phase 
uncertainty relation ANAcp > 1. The atomic electron/s during the two-photon excitation 
is/are acting as a device that measures the incident radiation phase to within an accuracy 
A(j) = u5t. Consequently the photon number would be uncertain by AA^ ~ ^, and this 
would shift the photoelectron spectrum by an amount ANhu), resulting in the suppression 
of AA^ lower-order peaks. This is a very well-known phenomenon in ATI (e.g. [18]) ; 
attempts to explain it have so far met with only partial success. Note that in the case of 
two electrons, they can act individually as such phase-measuring devices and no correlation 
is required ; as such, this feature is not restricted by the condition (24). Indeed, peak 
suppression (also known as quenching) has been observed at very high intensities (~ 10^^ 
W/cm2). 

Use of circularly polarised light lends a particularly interesting twist to this aspect. 
A circularly polarised photon can be regarded as a superposition of two linearly polarised 
components 90° apart in phase, and hence even if two circularly polarised photons arrive 
with a phase gap ivdt = ^ between them, the x-component of one photon and the y- 
component of the other can still be in exact phase correlation. This implies that the 

11 



coherence condition 5t <^ — can be relaxed to Jt ~ — , i.e. it would be satisfied at much 
smaller intensities. Thus at any given intensity, many more lower-order ATI peaks can 
be expected to be quenched by circularly polarised light than by linearly polarised light. 
This feature has indeed been observed (e.g. Bucksbaum et al [19], Bashkansky et al [20]). 
Interestingly, this also allows us to predict that with a circularly polarised laser, onset 
of DDI should occur in MPI experiments at lower intensities than with linearly polarised 
lasers. 

In conclusion, we find that our non-local model of QED gives us a correlated two- 
photon absorption which, unlike the conventional four-field correlation function, occurs in 
the first order of perturbation and is therefore linear in intensity. A single first-order matrix 
element [Ti, eqn.(17)] exactly reproduces the experimental finding of Georgiades et al [1], 
while in higher orders we obtain results that match very well with several experimental 
findings in MPI and ATI of atoms. From a detailed work-out of the non-local QED 
we obtain a non-thermodynamic arrow of time, and we also find that correlations in a 
quantum system must necessarily be of the EPR-type, a result that agrees with Bell's 
theorem. This non-local model QED can certainly be further improved upon by way of 
rigour and sophistication, but we believe that, basically, we have been able to demonstrate 
its usefulness. 
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APPENDIX A 

In this appendix we use relativistic notation, with units h = c = 1. Let A° and A 
be the scalar and vector potentials of an electromagnetic field. In field theory, the four- 
potential A^ = {A^,A) is usually introduced through imposition of the requirement of 
covariance under a U{1) gauge transformation (of the type e-*eA(x)-j^ ^^^ ^j^^g jg ^qhq via 
construction of the covariant derivative -D^ = dfj, + ieA^ {A^ = g^^jA^ , g^^, being the metric 
tensor) and replacement of the partial derivatives in the Lagrangian (e.g. of a complex scalar 
field (p) by these covariant derivatives. From this gauge-covariant Lagrangian D^(f)D^(f) we 
obtain the interaction Lagrangian —ej^A^, where j^ = i{(/)*d'^(j) — (^d^cj)*). (Of course, 
these A^'s need not be a dynamical field but can be a pure gauge field as well [21] ; however, 
this does not concern us here.) 

A significant point that is often lost sight of is that, in order to achieve gauge covari- 
ance, it is not a necessary condition that this interaction Lagrangian be local ; indeed, 
we can write down a non-local Lagrangian in a similar manner. Consider a non-local 
C/(l) gauge transformation of the kind exp(— ieA|x)(x'|). It is easily seen that a covariant 
derivative can be constructed as D^ = d^ + d^ + ieA^\x){x'\, (where dfj, = d/dx'^ and 
we have written A^\x){x'\ = g^^A'^\x){x'\), provided that A^ transforms under the gauge 
transformation as 

A^^ A^ + d^A\x){x'\ + df,A\x){x'\. (Al) 

It can be readily verified that this transformation property is consistent with the require- 
ment of gauge invariance of the non-local (electromagnetic) field tensor 

J''"'\x){x'\ = {d^ + d^')A''\x){x'\ - (a" + d'')A^\x){x'\ (A2) 

under this new gauge transformation, provided that 

[d^,d''] = 0. (A3) 

Ultimately, we have to replace A^ in the above-defined interaction Lagrangian by 
gtiuA'^\x){x'\. We see, therefore, that our non-local potential is consistent with the basic 
requirement of gauge covariance of the Lagrangian. Our eqn. (1) for the matrix element 
Mfi{t) follows directly from this non-local interaction. 

The nature of the correlation between x and x' in this formulation is laid down by 
the above commutation relation (A3), which indicates that x and x' must be completely 
independent, i.e. there is no restriction that their interval be time-like. We postpone a 
detailed discussion of this point till later ; in connection with eqn. (14) we shall show 
(see second paragraph after the eqn. in text) that energy conservation at the atomic level 
requires that this correlation must, of necessity, be of the EPR-type. 
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APPENDIX B 

We show in this appendix that only the retarded interaction can lead to energy con- 
servation. The time-integral in (14), with retarded correlation, is given by 

/CXD pt — p/c pOO pt — p/c 

^j-iEft-nut) ^^ ^-±iE^t'+n^t') ^^> ^ f^t)dt g{t')dt' (say) (Bl) 

-CXD J —OO J — OO J — OO 

{p = |r — r'l), while with the advanced correlation we have 

) /"CXD 

f{t)dt g{t')dt' (B2) 

D Jt-\-p/c 

The expression (Bl) can be written as a sum of two integrals : 

/O pt — p/c poo pt — p/c 

fit) dt / g{t') dt' + / fit) dt / g{t') dt' (Bl') 

-OO J — OO Jo J — OO 

while the expression (B2) can be written as 

/O /-OO ^oo poo 

f{t)dt g{t')dt'+ f{t)dt g{t')dt' (B2') 

-OO Jt+p/c Jo Jt+p/c 

Elementary considerations show that the first integral in (Bl') and the second integral in 
(B2') contribute nothing, while the other two integrals are (to within factors) given by 
S{Ef, Ei + 2huj). Energy conservation in the two-photon transition thus results, as it were, 
when we sum up our "samplings" of the phase of the matrix element over extended periods 
of time. 

We now come to the question of the significance of the limits of the integrals over 
dt. The matrix element is, so to say, created when we have photoexcited the atom ; let 
us agree to define that this act has been done at the instant t — 0. The crucial point to 
remember is that although the photoexcitation is not an irreversible process (the atom can 
certainly be de-excited again), our acts of exciting and subsequently observing the atom 
are irreversible in time. Therefore we can sample the phase (or whatever else we like to) 
of the excitation matrix element only at times t > 0. Thus the first integral in (B2'), 
although non-vanishing, is actually non-physical (it is just as meaningful as sampling the 
behavioural pattern of a baby before the baby is born), and we are left with only one 
non- vanishing integral in (Bl') 

poo pt — p/c 

/ e^^'^f'-^^'Ut e-^^''^''+^^''Ut' (B3) 

Jo J-oo 

Thus energy conservation will be satisfied if and only if temporal correlation is provided 
by the past, and not by the future ; we have derived an arrow of time in quantum physics 
as a necessary condition for energy conservation. 
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APPENDIX C 

Elementary analysis shows that 



poo poo 

TZ = Ri{r)r'^ dr I Rf{r)rdr, 

Jo Jo 

^1 = E ^^ E J KmA^^^nAr) dr j Y:Ar')Yi^^^ir') dr', 



n 1 



n=0 i^= — n ?Ti= — 1 

X [Y,Ur)y:Ar)yk^Ar)dr' 



The first angular integral gives If = k, while the second integral, which can be evaluated 
using the Gaunt formula, vanishes unless n = // ± 1 and n = /^ ± 1. Combining these two 
we get the selection rules // = h, ^i ± 2. 
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APPENDIX D 

The materials in this appendix are standard, and are included only for the sake of 
completeness. More details may be found in [9] (we have introduced some minor changes in 
notation). The Schrodinger equation for a two-electron atom in hyperspherical coordinates 
takes the form 

+ ^^] +V-E '^{R,a,f,r)=0 (Dl) 



2m V dR"^ i?2 
where the "grand angular momentum operator" 



d^_l /(/ + !) I'jl' + l) 
da^ 4 cos^ a sin^ a 



A^ = -^-7 + ^^S^ + ^^^^ (D2) 



Also, V = 7 + I — ^T, the total Coulomb interaction within the atom, and m is 

■ r r' |r— r'l ' ' 

the electron mass. The channel eigenfunctions $'s in (20) are obtained by configuration 
mixing in a basis of eigenfunctions of the grand angular momentum operator (D2). The 
eigenfunctions of this operator are 



■>nll'LM 



A{cos aY {sin af Pi""'^^ (cos 2a)y^ii' (r, r) (D3) 



where A is a normalisation constant, a = /' + ^, 6 = / + |, 

I I' 

yLii'= E E {Wmm'\WLM)YUf)Yi,^,if'), 

ra^ — l ra'^ — V 

Pi"''^(x) = i^2Fi(-n,a + 6 + n + l;a + l;i^) (D4) 

is the Jacobi polynomial, and n is the radial correlation quantum number, the same as rirc 
of Fano [9] . The eigenvalues of A^ are given hj \{\-\-l)% , where A = 2n + / + /'-!- |. Thus 
for each channel // in (20) we can write 

j 

where the running index j = (n, /, /', L, M) collectively. 

In hyperspherical coordinates the total Coulomb potential V is given by 

V=''^C{a,r,r'), C = -— — + .. . ./. . ,,,, (D5) 

n cos a sin a (1 — r.r sin iaj^/^ 

The hyperradial functions F^{R) in (20) are given by the coupled equations 

r^2 / ^2 



2m 



(-^ + ^^^) - e\ F,{R) + ^ $: V,,,F,,{R)S,l'Smm' 



(D6) 



Restricting ourselves to a single configuration in (20) and substituting eqn.(16) in (19), we 
obtain eqn. (21) for Tu. 
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